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ABSTRACT 

We use Keck/HIRES spectra of 37 optically bright quasars at z = 2-4 to study narrow absorption 
lines that are intrinsic to the quasars (intrinsic NALs, produced in gas that is physically associated 
with the quasar central engine). We identify 150 NAL systems, that contain 124 C iv, 12 N v, and 
50 Si IV doublets, of which 18 are associated systems (within 5,000 km s~-^ of the quasar redshift). 
We use partial coverage analysis to separate intrinsic NALs from NALs produced in cosmologically 
intervening structures. We find 39 candidate intrinsic systems, (28 reliable determinations and 11 that 
are possibly intrinsic). We estimate that 10-17% of C iv systems at blueshifts of 5,000-70,000 km s~^ 
relative to quasars are intrinsic. At least 32% of quasars contain one or more intrinsic C iv NALs. 
Considering N v and Si iv doublets showing partial coverage as well, at least 50% of quasars host 
intrinsic NALs. This result constrains the solid angle subtended by the absorbers to the background 
source(s). We identify two families of intrinsic NAL systems, those with strong N v absorption, and 
those with negligible absorption in N v, but with partial coverage in the C iv doublet. We discuss the 
idea that these two families represent different regions or conditions in accretion disk winds. Of the 
26 intrinsic C iv NAL systems, 13 have detectable low-ionization absorption lines at similar velocities, 
suggesting that these are two-phase structures in the wind rather than absorbers in the host galaxy. We 
also compare possible models for quasar outflows, including radiatively accelerated disk-driven winds, 
magnetocentrifugally accelerated winds, and pressure-driven winds, and we discuss ways of distinguishing 
between these models observationally. 

Subject headings: accretion - quasars: absorption lines — quasars: general — galaxies: active 



1. INTRODUCTION 

Quasars and active galactic nuclei (AGNs) are thought 
to be powered by accretion of matter onto a supermas- 
sive black hole. The accretion flow onto the black hole is 
thought to proceed via an equatorial accretion disk, which 
provides a mechanism for removing the angular momen- 
tum of the infalling matter. Outflows (winds) from such 
disks appear to be an inseparable part of this process. 
More specifically, hydromagnetic winds (e.g., Blandford 
& Payne 1982; Emmering, Blandford & Shlosman 1992; 
Konigl & Kartje 1994; Everett 2005) provide a potential 
mechanism for extracting angular momentum from the ac- 
creting material, thus allowing accretion to proceed. Ac- 
cretion disk winds are likely to be responsible for the broad 
absorption lines (BALs) observed in the spectra of a frac- 
tion of quasars (e.g., Murray et al. 1995; Arav, Li, & 
Begelman 1994, Proga, Stone, & Kallman 2000) and may 
even be the source of the broad emission lines that are the 
hallmark of all quasars and AGNs (e.g., Chiang & Mur- 
ray 1996; Murray & Chiang 1997). Understanding the 
quasar outflow mechanisms by testing and refining theo- 
retical models is therefore an essential part of our quest to 



understand the central engines of AGNs and quasars. 

Quasar outflows also have quite important consequences 
for cosmology and galaxy formation and evolution because 
they provide energy and momentum feedback to the in- 
terstellar medium (ISM) of the host galaxy and to the 
intergalactic medium (IGM). Simulations of galaxy evo- 
lution through mergers (Springel, Di Matteo & Hernquist 
2005) show that AGN feedback heats the ISM and inhibits 
star formation. Thus the colors of the resulting galax- 
ies evolve very quickly to the red, in agreement with the 
observed color distribution of nearby galaxies. In semi- 
analytic models of galaxy assembly (Granato et al. 2004; 
Scannapieco & Oh 2004) AGN feedback expells dense gas 
from the centers of the host galaxies, heating up the IGM 
(which it also enriches with metals). Because of its low 
density, the IGM does not cool efficiently and as a result 
it cannot fall back onto the galaxy and fuel star forma- 
tion. Thus feedback brings about an early termination 
of the assembly of the host galaxy. In this context it is 
quite important to know the fraction of quasars driving 
energetic outflows, as well as the energy content of the 
outflow. The latter quantity can be determined from the 
velocity, column density, and global covering factor of the 
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outflowing gas. All of these parameters are intimately con- 
nected to dynamical models of accretion disk winds and 
constraining them observationally can hardly be overem- 
phasized. Moreover, the outflowing material may provide 
a mechanism for enriching the IGM with metals. This idea 
is made more attractive by recent results that suggest high 
metallicities in quasar outflows (e.g., Hamann et al. 1997; 
Gabel et al. 2006; Petitjean et al. 1994; Tripp et al. 1996; 
D'Odorico et al. 2004). This question is closely connected 
to the issues above since the outflow rates and metallicities 
of quasar winds are at the heart of the matter. 

In the models cited above, the outflow is accelerated by 
either magnetoccntrifugal forces, or by radiation pressure 
in lines and continuum, or by a combination of the two 
processes. The flow originates deep inside the potential 
well of the black hole and can reach terminal speeds of 
order 10^ km s~^ or higher. The geometry of the flow dif- 
fers between a purely radiation-driven and a purely mag- 
netoccntrifugal wind. In the former case, (and especially 
at high luminosities) the flow is largely confined to low 
latitudes and effectively "hugs" the accretion disk. In ad- 
dition, the numerical simulations of Proga et al. (2000) 
show that the region above the fast, low-latitude stream 
develops transient filaments and streams which are denser 
than the ambient medium. In contrast, in the latter case 
the flow is nearly cylindrical (or U-shaped) and appears 
somewhat collimated at large radii from the black hole 
(sec, for example, the illustrations in Blandford & Payne 
1982 and Konigl & Kartje 1994). However, the density of 
the flow drops sharply once the gas has traveled a distance 
of a few launch radii, with the consequence that the dense 
part of the flow is concentrated near the equatorial plane 
(e.g., Everett 2005). A hybrid model, in which radiation 
pressure and magnetoccntrifugal forces are combined, also 
leads to a similar geometry (Everett 2005). One may ex- 
pect by analogy with radiation pressure-driven winds that 
the higher-latitude, lower-density regions of these winds 
would also contain transient, dense filaments or streams. 
Thus the apparent geometry and kinematics of the flows 
resulting from these two acceleration mechanisms are fairly 
similar. 

A different family of outflow models includes pressure- 
driven winds (e.g., Balsara & Krolik 1993; Krolik & Kriss 
1995, 2001; Chelouche & Netzer 2005). These are made up 
of gas that is photo-evaporated from the cool, dense torus 
that is invoked in AGN unification schemes. Thus, the flow 
originates not very deep inside the black hole potential well 
and the outflow speed is of order 10^ km s~^. The outflow 
is expected to have a multi-temperature and multi-density 
structure and indeed it would be clumpy. The primary dif- 
ference between this type of model and the accretion-disk 
wind models described above is the terminal speed of the 
flow. 

Blueshifted (or "P-Cygni") absorption lines provide a 
direct probe of quasar accretion disk winds. BALs (widths 
> 2,000 km s~^, by definition) have been the traditional 
means of studying such winds because they can be readily 
associated with them. This association is motivated by 
their large widths (often corresponding to velocity spreads 
of 30,000 km s~^) and their smooth profiles. However, 
intrinsic narrow absorption lines (intrinsic NALs; widths 
up to 500 km s~^; see review by Hamman & Sabra 2004) 



that are physically related to the quasars are an alternative 
and very useful means of studying such outfiows. Intrinsic 
(most often blueshifted) NALs are detected in a significant 
fraction (25-50%) of quasar spectra. These correspond to 
UV resonance transitions in highly ionized ions, such as 
C IV and N v (e.g., Foltz et al 1986, hereafter F86; Sargent, 
Steidel, & Boksenberg 1989; Anderson et al. 1987; Young, 
Sargent, & Boksenberg 1982; Steidel & Sargent 1991; Gan- 
guly et al. 2001, hereafter GOl; Vestergaard 2003). Quasar 
spectra also host an intermediate class of absorption lines, 
mini-BALs (widths between 500 and 2,000 km s~^, see 
Hamann & Sabra 2004), whose widths and smooth pro- 
files also suggest an origin in outflows. Such absorption 
lines are also fomid in Scyfcrt 1 galaxies (e.g., Crenshaw, 
Maran, & Mushotzky 1998; Crenshaw et al. 2004; Scott 
et al. 2004). Studies of intrinsic NALs and mini-BALs are 
more promising than studies of BALs for two important 
reasons: first these lines do not suffer from self-blending, 
and second, they are found in a wider variety of active 
galactic nuclei, and with greater ubiquity. In the context 
of the models described above, intrinsic NAL may probe 
the dense, cool filaments embedded in a hotter outflow. 
Thus they may give us access to a different portion, a dif- 
ferent phase, or a different line of sight through the outflow 
than the BALs. 

In spite of their potential importance, NALs have not 
received as much attention as BALs. Thus, we know sur- 
prisingly little about the relation between intrinsic NALs, 
mini-BALs, and BALs, and their connection, if any to the 
properties of the broad emission-lines. This is because 
it is difficult to distinguish intrinsic NALs from NALs 
that are not physically related to the quasars (interven- 
ing NALs), produced in intervening galaxies, intcrgalac- 
tic clouds. Milky Way gas or gas in the host galaxies of 
the quasars. It has been traditionally thought that NALs 
that fall within 5,000 km of quasar emission redshifts 
(termed associated absorption lines or AALs) are physi- 
cally associated with the quasars, because their frequency 
per unit velocity increases with decreasing velocity offset 
from the quasar (e.g., Weymann et al. 1979). AALs are 
common in all types of quasars, although the strongest 
ones appear preferentially in radio-loud quasar spectra 
(F86; Anderson et al. 1987). The statistical analysis of 
Richards et al. (1999) and Richards (2001) suggested that 
a fraction of 36% of the NALs with blueshifts from 5,000 
to 65,000 km s^^ from a quasar may also be physically 
associated with it. 

Over the past decade, with the advent of high-dispersion 
spectroscopy of faint objects, it has become possible to 
identify individual intrinsic NALs, based primarily on one 
or both of the following indicators: (a) the dilution of ab- 
sorption troughs by unocculted light (e.g., Hamann et al. 
1997a; Barlow & Sargent 1997; Ganguly et al. 1999, here- 
after G99; Ganguly et al. 2003; Misawa et al. 2003), and 
(b) time variability of line profiles (e.g., depth, equivalent 
width, and centroid), within a year in the absorber's rest 
frame (e.g.. Barlow & Sargent 1997; Wampler, Chugai, & 
Petitjean 1995; Hamann, Barlow, & Junkkarinen 1997b; 
Wise et al. 2004; Narayanan et al. 2004; Misawa et al. 
2005, hereafter M05). 

In this paper, we use the former of the above techniques 
to identiijy intrinsic C iv, Si iv, and N v NALs in the high- 
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resolution spectra of 37 quasars at Zem = 2-4. These 
quasars were originally selected without regard to NAL 
properties, though BAL quasars were avoided, and there 
was a preference for optically bright quasars. These data 
allow us to construct a large, and relatively unbiased, sam- 
ple of intrinsic NALs. Using this sample, we investigate the 
NAL properties statistically and we compare them with 
the properties of the quasars that host them. Since our 
spectra typically cover several transitions of the same sys- 
tem, we are also able to probe the ionization state of the 
absorber. This large sample represents a major advance 
over most previous efforts, which either dealt with small 
samples of intrinsic NALs or employed low- resolution spec- 
tra in statistical studies. 

In we describe the properties of the quasar sample 
and briefly summarize the observations. In [J3]we describe 
our methodology for identifying intrinsic NALs and evalu- 
ating the reliability of this determination. Our results are 
presented in Sj4l their implications in the context of 
models for quasar outflows arc discussed in fj5l Our con- 
clusions are summarized in SJH] In this paper, we use a cos- 
mology with Ho—75 km s^^Mpc^^, 17^=0.3, and ^1a=Q.7. 
In the Appendix we present all the detailed information on 
the NALs we have detected. Included in this Appendix are 
(a) a table of NAL properties derived from fits to their pro- 
files, (b) diagnostic plots of UV resonance doublets, which 
form the basis of our intrinsic NAL identification method, 
and (c) comparison plots of profiles of different transitions 
from the same intrinsic NAL system, (d) detailed notes on 
individual intrinsic NAL systems. 

2. QUASAR SAMPLE AND OBSERVATIONS 

The quasars in our sample were originally selected and 
observed in a survey aimed at measuring the deuterium-to- 
hydrogen abundance ratio (D/H) in the Lya forest. The 
typical value of D/H is so small, 2-4x 10~^ (O'Meara et al. 
2001 and references therein), that we can detect only D i 
lines corresponding to H i lines with large column densities, 
log(iVH i/cm~^) > 16.5. Therefore the survey included 40 
quasars, in which either damped Lya (DLA) systems or 
strong Lyman limit systems (LLSs) were detected. The 
observations were carried out with Keck/HIRES through 
a 1'.'14 slit resulting in a velocity resolution of ~ 8 km s^^ 
(FWHM). The spectra were extracted by the automated 
program, MAKEE, written by Tom Barlow. In this paper 
we use the spectra of 37 of these quasars, listed in Table[Tl 
which cover the rest-frame wavelength range between the 
Lya and C IV lines. This target selection method does 
not directly bias our sample with respect to the properties 
of any intrinsic absorption-line systems that these spectra 
may contain. However, an indirect bias could result since 
the sample contains only quasars bright enough to allow 
high S/N, high resolution observations. Although, opti- 
cal brightness does not appear to be the most significant 
factor in determining whether a quasar hosts a NAL, it is 
likely to have some effect (GOl). Also, there is a some- 
what enhanced probability of finding associated NALs in 
quasars that host BALs (GOl), thus the avoidance of BAL 
quasars in our sample may bias against intrinsic NALs. 
We will discuss these possible biases, and comparisons to 
other samples selected by different criteria, in § 14.11 

Quasar emission redshifts were obtained from a variety 



of sources in the literature and they are based primar- 
ily on measurements of the peaks of their strong, broad 
UV emission lines, namely Lya, Si iv and C iv. It is 
well known that the redshift determined from these par- 
ticular lines is systematically different from the redshift 
of the low- ionization lines (i.e., Mg ii and the Balmer 
lines) and the systemic redshift, as given by the narrow, 
forbidden lines (see, for example, Corbin 1990; Tytler & 
Fan 1992; Brotherton et al. 1994; Sulentic et al. 1995; 
Marziani et al. 1996). Even though redshift differences 
can reach ~ 4, 000 km s^^, in most cases they are within 
±1, 000 km s^^. In radio-quiet quasars, there is a system- 
atic tendency for the C iv line to be blueshifted relative 
to the systemic redshift. Most relevant to this work are 
the results of Tytler & Fan (1992) who study the redshift 
discrepancies of the UV lines in quasars of comparable red- 
shifts and luminosities to those of our sample. They find 
a mean blueshift of the UV lines relative to the systemic 
redshift of 260 km s~^ and that 90% of the blueshifts are 
between and 650 km s^^ (see their Figure 13). 

We also list in Table [T] a measure of the radio- loudness of 
the quasars in our sample, the ratio of the flux densities at 
5 GHz and 4400 A, i.e., 7^ = /^(5 GHz)//j.(4400 A) (fol- 
lowing Kellermann et al. 1989, 1994). We adopt 7^ > 23 
as the criterion for radio-loudness (instead of the TZ > 10 
adopted by Kellermann et al.) in order to separate cleanly 
the low-luminosity radio sources in our sample from the 
high-luminosity ones, which have TZ > 252. We derived 
fi, (4400 A) from the V or R magnitude using the following 
equations (Schmidt & Green 1983; Oke & Schild 1970), 

m,/ = -2.51og/i.(5500 A) - 48.60 (1) 

rrifl.^ -2.5 log/^(6600 A) -48.82 (2) 

and assuming an optical spectral index of ao — 0.44 (where 
oc I/-"; see Vanden Berk et al. 2001). For Q0241-0146 
(zem = 4.040) and Q1055-h4611 {z^m = 4.118), the ob- 
served optical fluxes obtained from the V magnitude are 
underestimated because of contamination by the Lya for- 
est at Zabs~ 3.5. Therefore, we boosted these fluxes by 
dividing by the normalized transmission in the Lya forest 
at redshift z, 

T = exp [-0.0037 (1 -I- zf-'^^] (3) 
(Press, Rybicki, & Schneider 1993). The value of 
/^(5 GHz) was obtained from measurements of the flux 
density at various radio frequencies, assuming a radio spec- 
tral index of a^ = 0.7. Thus, our 37 quasars are separated 
into 12 radio-loud quasars and 25 radio-quiet quasars. Un- 
fortunately, the number of radio-loud quasars that cover 
the necessary wavelength regions for our analysis is too 
small to allow a useful statistical study of differences in 
NAL properties between subclasses. 

The properties of the quasars in our sample are sum- 
marized in Table [Hand Figure [T] Columns (1) and (2) of 
Table [1] give the quasar name and emission redshift, and 
columns (3) and (4) the V and R magnitudes. The optical 
flux density at 4400 A in the rest frame is given in col- 
umn (5). Column (6) lists the observed radio flux at the 
frequency of column (7), if a radio source was detected 
within 10" of the optical source (e.g., Kellermann 1989). 
The derived radio flux at 5 GHz is listed in column (8). 
Column (9) gives the radio-loudness parameter, TZ, based 
on which the quasars are labeled as L (radio-loud) or Q 
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(radio-quiet) in column (10). The distributions of optical 
(at 4400 A) and radio (at 5 GHz) luminosities are plotted 
in Figure [11 

3. SAMPLE OF NAL SYSTEMS 

To construct our sample of NAL systems for statisti- 
cal analysis, we first examine the 37 quasar spectra, and 
mark all absorption features that are detected at a confi- 
dence level greater than 5a, i.e., [1 — Rc,obs]/cy{Rc,obs) > 5 
{Rc,obs is the observed residual intensity at the center of a 
line in the normalized spectrum and cr{Rc^obs) is its uncer- 
tainty). Next, we identify N v, C iv, and Si iv doublets in 
the following regions around the corresponding emission 
line: 

N V absorption doublets: from —5,000 km s"""^ to 
km around the N v emission line in the spec- 
tra that include this line0 If the spectrum extends 
redward of the N v emission line, we also search 
for absorption doublets up to H-10,000 km s^^ to 
the red of the line. The velocity range is relatively 
narrow for this transition because of contamination 
from the Lya forest, which is rather serious at the 
redshift of our target quasars. 

C IV absorption doublets: from -70,000 kms~^ to 
km s~^ around the C iv emission line in the spec- 
tra that include this line. If the spectrum extends 
redward of the C iv emission line, we also search 
for absorption doublets up to -f 10, 000 km s^^ to 
the red of the line. 

Si IV absorption doublets: from —40, 000 km to 
km s^^ around the Si iv emission line in the spec- 
tra that include this line. If the spectrum extends 
redward of the Si iv emission line, we also search 
for absorption doublets up to -|-10,000 km s^^ to 
the red of the line. 



In Table [U we list the red limit of the velocity window 
that we searched for each quasar. 

Absorption troughs that are separated by non-absorbed 
regions are considered to be separate lines, even if they are 
very close to each other. The equivalent width is measured 
for each separate line by integrating across the absorption 
profile. In the case of doublets, we represent the equiv- 
alent width by that of the stronger (blueward) member. 
In total, 261 C iv, 13 N v, and 92 Si iv doublets are 
identified in this manner. To facilitate detailed studies of 
the systems, we also searched for single metal lines at the 
same redshifts as the doublet lines even if these were in 
the Lya forest. For the statistical analysis, we construct a 
complete sample that contains only doublet lines whose 
stronger members would be detected even in our most 
noisy quasar spectrum. We convert the absorption-feature 
detection criterion given above, [1 — Rc.obs]/ (^(Rc,obs) > 5, 
to an equivalent width limit using equation (1) of Mis- 
awa et al. (2003) and we select absorption features above 
a certain rest-frame equivalent width. The rest-frame 
equivalent width limits are set by the spectrum of the 
quasar Q1330-h0108 as follows: W„,„ (C iv) = 0.056 A 

^ The velocity offset is defined as negative for NALs that are blueshifted from the quasar 
As we noted in JS] a "line" is separated from its neighbors by non-absorbed regions. 



at A « 5800 A, Wm,n (N v) = 0.038 A at A « 5700 A, 
and W^^n (Si iv) = 0.054 A at A « 5800 A. These limits 
apply to the stronger member of each doublet. The "ho- 
mogeneous" NAL sample defined by these limits contains 
138 C IV, 12 N V, and 56 Si iv doublets. 

In order to evaluate the column densities [N in cm~^) 
and Doppler parameters (6 in km s~^) of NALs, we need to 
deblend the absorption profiles into narrower components. 
Thus, we used the software package minfit (Churchill 
& Vogt 2001) to fit absorption lines with Voigt profiles. 
In the fitting process, we consider the coverage fraction 
(defined and discussed in detail in ij3.11 below) as a free 
parameter, as well as the redshift, column density, and 
Doppler parameter, and we convolve the model with the 
instrumental profile before comparing with the data. Kine- 
matic components are dropped if a model with fewer com- 
ponents provides an equally acceptable fit to the data. 
With this procedure, the NALs of the homogeneous sam- 
ple are deblended into 483 C iv, 41 N v, and 182 Si iv 
components. The observed profiles of C iv, N v, and Si iv 
NALs, best-fitting models, and la error spectra, are shown 
on a velocity scale relative to the flux weighted line center 
in the figures of the Appendix. Almost all the NALs are 
deblended into multiple components. 

To refine our sample for statistical analysis, we combine 
NALs that lie within 200 km s~^ of each other into a sin- 
gle system (a so-called "Poisson system"). This method 
is based on the assumption that clustered lines are not 
physically independent (e.g., Sargent, Boksenberg, & Stei- 
del 1988). We chose this clustering velocity after con- 
structing the distribution of velocity separations between 
C IV NALs, which we show in Figure [21 We note that 
our adopted clustering velocity differs from the value of 
1,000 km s~^ adopted in previous papers (e.g., Sargent et 
al. 1988; Steidel. 1990; Misawa et al. 2002). We obtain 
150 Poisson systems in this manner, of which 124, 12, and 
50 systems contain C iv, N v, and Si iv NALs, respectively. 
All NALs are also classified as AALs and non-AALs. We 
found 18 associated Poisson systems, of which 9, 12, and 
4 systems contain C iv, N v, and Si iv NALs, respectively. 

Hereafter, we will use the terms "(Poisson) system" for 
a group of NALs within 200 km s~^ each other, "line" 
or "doublet" for an individual NAL (that may lie within 
a Poisson system), and "component" for a narrow kine- 
matic component described by a single Voigt profile (as 
deblended by minfit In Table [H we list the total num- 
bers of systems, lines, and components that we have found 
for each of the three strong absorption lines of interest, 
namely C iv. Si iv, and N v. 

We compute the flux- weighted center of a line (or a Pois- 
son system) as the first moment of the line profile, i.e., 

EA. [l-R{K)] 



^obs — (A) — 



E [1 - 



(4) 



and the line dispersion, s, via 



where (A^) 

by 



= v/(A2> - (A)2 , (5) 
is the second moment of the line profile, given 



(A^ 



EAf [l-i^(A.) 

E[i-^(A.)] 



(6) 
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The above sums are evaluated over the entire width (fuU 
width at zero depth) of a line. The quantity R{Xi) is the 
intensity at pixel i in the normalized spectrum and Xi is 
the wavelength of that pixel. Using the above quantities, 
we obtain the flux-weighted line width from 



cr(w) = \/2 s„ = a/2 ^-^^ c (km s 



(7) 



Using the flux-weighted wavelength of a line we compute 
its flux- weighted redshift, and then its velocity offset by 
the relativistic Doppler formula, i.e., 

_ Vshift _ (1 + Zem)^ - (1 + Zabs)'^ 



(3: 



(8) 



C {I ^ Zemf + {I + Zabsf ' 

where Zem and Zabs are the emission redshift of the quasar 
and the absorption redshift of the line, respectively. The 
wavelengths, redshifts, velocity offsets, and flux-weighted 
widths of detected lines, computed as described above, are 
given in the Appendix, where we also include NALs that 
do not meet the rest-frame equivalent width criterion for 
our homogeneous sample (i.e., they have Wrest < Wmin)- 

3.1. Partial Coverage Analysis 

We identified intrinsic NAL candidates by looking for 
absorption troughs that were diluted by unocculted light 
from the background source. The optical depth ratio of 
doublet lines, such as C iv. Si iv, or, N v, sometimes de- 
viates from the value of 2:1 expected from atomic physics. 
This discrepancy can be explained if the absorber covers 
the background flux source only partially, and the unab- 
sorbed flux changes the relative depth of the lines (e.g., 
Wampler et al. 1995; Barlow & Sargent 1997; Hamann et 
al. 1997a; G99). However, there are other possible expla- 
nations, such as local emission by the absorbers (Wampler, 
Chugai, & Petitjean 1995) and scattering of background 
photons into our line of sight (G99). We assume the sim- 
plest model, in which the absorber has a constant opti- 
cal depth across the projected area (i.e., a homogeneous 
model). Then, the observed intensity as a function of ve- 
locity from the line center in the normalized spectrum is 
given by 

R{v)^[l-Cf{v) + Cf{v)e-^^^^] , (9) 

where Cf{v) is the fraction of background light occulted by 
the absorber (hereafter, the "coverage fraction"), and t{v) 
is the optical depth at velocity w. If a NAL system consists 
of many kinematic components (the most common case), 
these are combined by multiplying the individual Voigt 
profiles after making the appropriate adjustments for par- 
tial coverage, following equation Thus we describe the 
final normalized residual intensity by the product 

Riv) = n [i - CfMv) + CfMv) ' (10) 

k 

where the index k labels different kinematic components. 
This assumption is quite safe if the components do not 
overlap significantly in velocity, and is still a good approx- 
imation if one component dominates at each wavelength. 

If we measure the optical depth ratio of doublet lines 
with oscillator strength values of 2:1 (e.g., C iv, N v, and 
Si iv) by fitting their profiles, we can evaluate the coverage 
fraction as a function of velocity across the line profile as 



where Rr and Rb are the continuum normalized intensi- 
ties of the weaker (redder) and stronger (bluer) members 
of the doublet (see, for example. Barlow & Sargent 1997). 
G99 refined this technique by considering two background 
sources: the continuum source and the broad emission line 
region (hereafter BELR) . In this composite picture the to- 
tal coverage fraction can be expressed as a weighted aver- 
age of the coverage fractions of the two regions, namely. 



Cc{v) + w{v) Cbelb.{v) 

^ iTW) ■ 



(12) 



Cfiv) 



[l~Rr{v)Y 



l + Rb{v)^2Rr{v)' 



(11) 



where Cc(w) and Cbelr{v) are the coverage fractions of 
the continuum source and the broad emission line region, 
and w{v) is the ratio of the broad line flux to the contin- 
uum flux at a given pixel in the absorption-line profile. 

It is quite possible for the absorber to have different op- 
tical depths along different paths (i.e., an inhomogeneous 
model; deKool, Korista, & Arav 2002). Nonetheless, by 
exploring this inhomogeneous model, Sabra & Hamann 
(2005) found that the average optical depths derived from 
homogeneous and inhomogeneous models are consistent 
with each other within a factor of <l-5, unless a small 
fraction of the projected area of the absorber has a very 
large optical depth. Moreover, in most cases, one can 
find acceptable fits to the absorption profile using a simple 
combination of just a few Voigt components. Therefore, 
we adopt the homogeneous optical depth model, as de- 
scribed by equation (jlip. with the understanding that C/ 
describes the fraction of all background photons that pass 
through the absorber. We consider the composite picture 
described by equation (|12p . only if it critically affects our 
conclusions. 

We use two methods to evaluate C/; (i) the pixel-by- 
pixel method (e.g., G99), in which we apply equation pT|) 
to each pixel in the normalized spectrum, and (ii) the fit- 
ting method (e.g., Ganguly et al. 2003) in which we fit the 
absorption profiles using MINFIT, treating C/ as well as 
logiV, b, and Zabs as a free parameter. In the former case 
we obtain a value of C/ for each pixel in the line profile, 
while in the latter case we obtain a value of C/ for each 
kinematic (Voigt) component. Our derived C/ values are 
subject to a number of additional uncertainties resulting 
from line blending, uncertainties in the continuum fit, the 
convolution of the spectrum with the line spread function 
(LSF) of the spectrograph (G99; only apphes to the pixel- 
by-pixel method), and Poisson errors in regions with weak 
lines (only applies to the fitting method). The method- 
ology for assessing these uncertainties was developed in 
M05 with the help of extensive simulations. We do not 
use spectral regions that are strongly affected by these un- 
certainties, and we consider these sources of error when 
evaluating whether C/ deviates from unity, as described 
below. The C/ values evaluated by both of the above 
methods are overplotted on the observed spectrum of each 
doublet shown in Figure [T^] of the Appendix. 

The MINFIT routine sometimes gives unphysical cover- 
age fractions such as C/ < or C/ > 1. Tests by M05 
showed that coverage fractions produced by minfit are 
very sensitive to continuum level errors, especially for very 
weak components whose C/ values are close to 1, which 
suggests a cause for these unphysical results. In such 
cases, other fit parameters, such as the column density and 
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Doppler parameter, have no physical meaning. Therefore 
we re-evaluate them assuming C/ = 1 (with no error bar) 
for those components, and by refitting the rest of the line. 

Based on the results of the coverage fraction analysis, we 
separate all NALs into three classes based on the reliabil- 
ity of the conclusion that they are intrinsic: classes A and 
B respectively contain "reliable" and "possible" partially 
covered NALs (i.e., intrinsic NAL candidates), and class 
C consists of NALs with no evidence for partial coverage 
(i.e., intervening or unclassified NALs). We may regard 
classes A-C as including progressively smaller fractions of 
intrinsic NALs and we will treat them as such in our sub- 
sequent statistical analysis. 

We classified four mini-BALs as class A systems. We 
also classified three NALs into class B, although they 
did not exhibit partial coverage, because they show line- 
locking, which is often seen in intrinsic NALs at Zabs ~ 
Zem ■ In the list below and in the flow chart of Figure [31 
we give the criteria and method that we used to classify 
NALs into reliability classes. 

Class A: Reliable intrinsic NAL candidates: 

Al: smooth and broad (self-blended) line profile 
(i.e., mini-BAL). 

A2: MINFIT gives Cf+ 3a{Cf) < 1 for at least one 
component. 

Class B: Possible intrinsic NAL candidates: 

Bl: Both MINFIT and the pixel- by- pixel method 
give Cf+ <j{Cf) < 1 at the center of a com- 
ponent for at least one kinematic component. 

B2: line-locked. 

Class C: Lines without partial coverage and unclassifi.- 
able lines: 

CI: MINFIT gives C/ -I- <j{Cf) > 1 for aU compo- 
nents. 

C2: No components can be used for classification 
because they suffer from systematic errors as 
described in the next paragraph. 

C3: The system is not acceptable for classification 
because of problems due to model fit, contin- 
uum fit, or critical data defect. 

During line classification, we ignore components that 
satisfy the following criteria. These were established by 
following the same procedure to fit strong Mg ii doublets 
identified in the same spectra. For the strong Mg ii lines 
we would expect full coverage for all components, but we 
find C/ to be inconsistent with unity if these criteria apply: 

1. Weak components at the edge of a system. 

2. Weak components between much stronger compo- 
nents. 

3. Components in a heavily blended region (except for 
regions with extremely high S/N). 



4. Components for which we obtained unphysical val- 
ues of C/ in the first fitting trial. As we noted 
above, we set C/ = 1 for such components and re- 
peat the fit for the other components in the system. 

If at least one component of a NAL satisfies one of the 
criteria for class A or B, we would classify the NAL as 
an intrinsic NAL candidate, even if the other components 
are consistent with full coverage. This applies even for a 
Poisson system for which only one component in one of 
the separate lines satisfies the class A or B criteria. Also, 
if a Poisson system has more than one doublet transition 
detected and one or more of the transitions (C iv, N v, or 
Si iv) belongs in class A, it is considered a class A intrinsic 
system. Similarly, such a system will be considered a class 
B system if one or more transitions fall in class B (and 
none is class A). 

Placing a Poisson system in class A or B if only one 
component satisfies the classification criteria involves the 
following caveat: the class A or B NAL may be superposed 
on an intervening system by chance. We have thus eval- 
uated the probability of this chance superposition, using 
the information in Tabled For example, we have searched 
for associated C iv NALs in a velocity range Sf3 = 0.17 
and found 3 class A-l-B and 6 class C NALs. Since each 
NAL typically spans a velocity window of 200 km s~^ or 
less, the class C NALs cover 2.4% of this velocity window. 
Therefore, the probability that any one of the 3 class A-l-B 
NALs would be blended with a class C NAL (i.e., the ex- 
pected number of blends) is 0.07. Similarly, we estimate 
the number of blends to be 0.25 in the non-associated C iv 
velocity range, 0.03 in the N v velocity range (associated 
only), and 0.07 in the Si iv velocity range (non-associated 
only). Thus, we conclude that the danger of a blend is 
negligible. 

In our final sample of 150 Poisson NAL systems, 28 fall 
into class A, 11 into class B, and 111 into class C. The 
number of systems, lines, and components are summarized 
in Table O In the Appendix, we list the detailed classi- 
fication results, along with the subcategory numbers that 
give the basis for the classification (e.g., A2, C3, ...). In the 
Appendix, we also give detailed comments on individual 
NAL systems. In Tabled we give a census of the resulting 
Poisson systems, NALs and kinematic components. 

4. RESULTS 

Using the complete NAL sample constructed in the pre- 
vious section, we first investigate the density of NALs per 
unit redshift interval and per unit velocity interval (i.e., the 
velocity offset distribution, dN/dz and dN/d(3 = c dN/dv, 
respectively). Next, we examine the relative numbers of 
intrinsic NALs in different transitions. This is of inter- 
est because these relative numbers could be indicative of 
the ionization structure of the absorbers and also of their 
locations relative to the continuum source. Our analysis 
provides only lower limits on the intrinsic NAL fractions, 
because intrinsic absorber need not always to show partial 
coverage, and because some of our sample quasar spectra 
cover only very small redshift windows for NAL detection. 
We also examine the distributions of coverage fractions 
and line profile widths, and consider the ionization con- 
ditions of the intrinsic NALs in our sample. Finally, we 
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briefly compare the properties of quasars and the proper- 
ties of intrinsic NALs. 

4.1. Velocity Offset and Equivalent Width Distribution 

In evaluating dN/dz, we exclude segments of spectra in 
echelle order gaps, which affect about 8.8%, 16.0%, and 
23.0% of the spectral regions at ~5,500 A, ~6,000 A, 
and ^^6,500 A, respectively. We implicitly assume that 
all NALs arise in outflows from the quasars and we use 
the velocity offset computed by equation ([8]). In Table [3l 
we summarize our derived values of dN/ dz and dN/dP for 
different transitions and for different velocity offsets rela- 
tive to the quasar, and we also break down the results by 
reliability class. The velocity offset distributions of C iv, 
N V, and Si iv NALs are presented in Figure 21 broken 
down according to reliability class. It is noteworthy that 
intrinsic NAL candidates (i.e., class A and B) are found 
not only among AALs but also among non-AAL systems. 

In previous studies using low resolution spectra, the ve- 
locity offset distributions of C iv NALs are found to be 
almost uniform up to Vshift ^ —70, 000 km s^^ (the max- 
imum velocity at which the C iv doublet can be detected 
without blending with the Lya forest). However, a sig- 
nificant excess of AALs was found within 5,000 km s~^ 
of radio-loud and steep-spectrum quasars (Young, Sar- 
gent, & Boksenberg 1982; Richards et al. 1999, 2001; 
F86). Weymann et al. (1979) presented evidence that 
the distribution of intrinsic C iv NALs extends up to 

Vshift 18,000 km s-i. Richards et al. (1999) found 

an excess of C iv NALs with v shift ^ — 65, 000 km s~^ in 
optically luminous quasars, in radio-quiet over radio-loud 
quasars, and in flat-radio spectrum quasars over steep- 
radio spectrum quasars. Our results show no remarkable 
excess of NALs in radio-loud quasars, nor any strong ex- 
cess of NALs near the quasars in either the radio-loud or 
quiet subsamples. This is not surprising, however, since 
our subsamples of radio-loud and radio-quiet quasars are 
relatively small and heterogeneous. 

Considering all C iv NALs (from classes A, B, and C) 
together, we do find (see Table [3]) that, for C iv NALs 
the values of dN/dz and dN/df] in associated regions are 
about twice as large as those in non-associated regions, 
in general agreement with the results of Weymann et al 
(1979). In Figure O we show the velocity offset distribu- 
tion of all C IV NALs from our survey compared to that 
of class A-l-B NALs. 

The distribution of rest-frame equivalent widths of the 
C IV Poisson systems detected in our survey is shown in 
Figure [HI This distribution rises sharply towards our de- 
tection limit of 0.056 A. We note that the vast majority 
of our NALs are considerably weaker than those included 
in previous surveys at low spectral resolution. For exam- 
ple, the survey of Weymann et al. (1979) has a rest frame 
equivalent width limit of about 0.6 A (higher than 87% of 
the NALs in our C iv sample) while the surveys of Young 
et al. (1982) and F86 have a detection hmit of 0.3 A 
(higher than 70% of the NALs in our sample). The more 
recent Vestergaard (2003) study was also conducted using 
low spectral resolution, so it focused on C iv absorption 
lines stronger than ^0.3 A. 

A comparison of the equivalent width distributions be- 
tween our study and these previous low resolution stud- 



ies is a useful indication of possible biases due to quasar 
properties. However, such a comparison must focus only 
on our stronger systems. Also, we must group together 
different lines that fall within 200 km s"^ of each other 
in order to compare with systems found at low resolution 
(as we did in Figure [S]) . First we focus on systems with 
Wr{C iv) > 1.5 A since we are concerned about possible 
biases due to rejecting BAL quasars from our sample, and 
since Vestergaard (2003) targeted the strongest systems as 
likely intrinsic candidates. If we take into account Vester- 
gaard's different method of measuring the equivalent width 
of C IV, we find that her sample includes 6/114 quasars 
that have systems with Wr{G iv) > 1.5 A. However, the 
Vestergaard (2003) sample has approximately equal num- 
bers of radio-loud and radio-quiet quasars, unlike ours, in 
which only 1/3 of the quasars are radio- loud. If we con- 
sider only her radio-quiet sample, only 1/48 quasars has a 
C IV NAL with Wr{C iv) > 1.5 A. In her radio-loud sam- 
ple, 5/66 quasars have very strong C iv NALs. Thus we 
would predict that in our sample of 37 quasars (with only 

12 radio-loud) we would find 1.4 very strong C iv NALs. 
Though we did not find any, our number is certainly con- 
sistent with this estimate, and we do find a system with 
Wr{C iv) ~1.0 a, just below that hmit. We note that 
due to the nature of Vestergaard's study, her sample has a 
larger than representative fraction of radio-loud quasars. 
Our sample also has a somewhat larger fraction of radio- 
loud quasars (37%) compared to the general population 
(~10%), which would lead to a bias towards detecting very 
strong NALs compared to the general population. 

We next consider somewhat weaker NALs, down to 
Wr{C iv) = 0.3 A, to facilitate another comparison with 
the work of Vestergaard (2003). She found that 27% of 
the quasars in her sample had C iv absorption for which 
the sum of the equivalent widths of the two members of 
the doublet were > 0.5 A. In our sample, we only have 

13 quasars for which the most of the C iv associated re- 
gion is covered. Of those 13 quasars, 9 have Wr{C iv) 
> 0.056 A NALs, but only 1(8%) has Wr{C iv) > 0.3 A. 
This is somewhat smaller than we would predict (3.5/13 
= 27%) based on Vestergaard's results, however we again 
note the differences in the radio properties between the 
samples. Furthermore, we do not have full coverage of the 
associated region in all of the 13 quasars, and only 3 of 
the 13 are radio-loud. It is also worth noting that typi- 
cally the Vestergaard (2003) quasars are 1-2 magnitudes 
less luminous than the quasars in our sample. 

4.2. Line Width Distribution 

Historically, it has been very difficult to separate in- 
trinsic NALs from intervening NALs without high resolu- 
tion spectra because of the similarities of their line pro- 
files. However, if intrinsic NALs arise from absorbers 
similar to those producing BAL and mini-BAL systems 
(widths >, 2,000 km s~^), they could have larger widths 
compared to those of intervening NALs. If a tendency for 
large widths is found among intrinsic NALs, this may be 
used as an additional indicator of their nature. In Fig- 
ure [71 we compare the distribution of the fiux-weighted 
line width [ct(u); see eqn. (I?])] of intrinsic NALs (classes 
A and B) to that of intervening/unclassified NALs (class 
C). We do not see any substantial differences between the 
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distributions that would allow us to distinguish between 
intervening and intrinsic NALs. The average value of a{v) 
for intrinsic NALs (53.4 km s~^) is slightly larger than 
that for intervening NALs (44.5 km s~^). A Kolmogorov- 
Smirnov test gives a probability of 12% that the two distri- 
butions have been drawn by chance from the same parent 
population. 

4.3. Relative Numbers of Intrinsic NALs in Different 
Transitions 

The numbers of C iv, N v, and Si iv NALs in each of the 
coverage fraction classes are listed in Table We also list 
in Table [3] the values of dN/ dz and dN/d(3 for each type 
of system. Based on these values, reliable (class A) intrin- 
sic NALs make up 11% of all C iv NALs, 75% of all N v 
NALs, and 14% of all Si iv NALs. These fractions increase 
to 19%, 75%, and 18% after adding possible intrinsic NALs 
(class B). Concentrating on different transitions in turn, 
we note the following: 

N V. — The above fractions of N v NALs that are intrin- 
sic refer only to the AAL velocity range {v shift < 
—5, 000 km s^^) because contamination by the Lya 
forest prevented us from searching for NALs at 
larger velocity offsets. Nonetheless, it is worth not- 
ing that 3/4 of these associated N v NALs are in- 
trinsic. The density of intrinsic NALs in the asso- 
ciated region is thus quite high {dN/dz = 3.4 and 
dN/df3 = 13.7). 

Civ. — The fraction of intrinsic C iv NALs in asso- 
ciated regions (33%) is slightly higher than the 
fraction we found in non-associated regions (10%- 
17%). If the fraction of intrinsic NALs at Vshift > 
—5,000 km was actually the same as the frac- 
tion at Vshift < —5,000 km s"^, the probability 
that 3 or more out of 9 C iv AALs are intrinsic 
to the quasar (based on the binomial distribution) 
is 5.1% . If we consider only the radio-quiet sub- 
sample, the probability is even smaller, 1.2%. In 
the radio-loud sample, we found no intrinsic C iv 
AALs, but this does not mean that there is an ac- 
tual deficiency of intrinsic C iv AALs, because our 
12 spectra of radio- loud quasars have only a small 
redshift path for detecting C iv AALs {Sz = 0.64). 

Si IV. — We found no intrinsic Si iv AALs in the redshift 
window of Sz = 1.26 that we searched. 

4.4. Fraction of Quasars Showing Intrinsic NALs 

To investigate the geometry of absorbing gas around the 
quasars we count the number of intrinsic NALs per quasar, 
and we evaluate the fraction of quasars hosting intrinsic 
NALs. These values constrain the global covering factor 
(r2/47rj3 of the continuum source and BELR by the ab- 
sorber and the distribution of absorbing gas around the 

* We considered also subsamplcs of radio-loud and radio-quiet quasars but found no significant differences from the total sample. 
^ Here it is important to make the distinction between the terms "coverage fraction" and "global covering factor" . The former term (defined 
in refers to the fraction of the projected area of the background source occulted by an absorbing parcel of gas or, equivalently, to the 

fraction of background photons that pass through the absorber. The latter term refers to the total solid angle subtended by the ensemble of 
absorbers around the source. 

A parcel of gas in the vicinity of the continuum source could be large enough to cover the source completely. In such a case our partial 
coverage method will yield Cf = I, even though the absorption lines are really intrinsic. 



quasar central engine. At z < 1, GOl found that radio- 
loud, flat-spectrum quasars with compact radio morpholo- 
gies (i.e., with a face-on accretion disk) lack AALs down 
to a limiting equivalent width of Wrest = 0.34 A. This 
suggests that the presence of intrinsic NALs depends on 
the inclination of the line of sight to the symmetry axis of 
the quasar central engine. 

In Table [3] we list the fraction of NALs that are intrinsic 
and the fraction of quasars that have intrinsic NALs. We 
count only quasars that have velocity windows for NAL 
detection of S/3 > 0.01 {Sv > 3,000 km s~^) in the tran- 
sition of interest. In Figure [51 we plot the distribution 
of the number of intrinsic NALs per quasar. The upper 
panel shows the distribution of only class A NALs, while 
the lower panel includes both class A and B NALs. 

About 24% of quasars (9 out of 37) have at least one re- 
liable (class A) intrinsic C iv NAL. Taking possible (class 
B) NALs into consideration, 32% (12 out of 37) of quasars 
have at least one intrinsic C iv NAL. The corresponding 
fractions for N v and Si iv (19-24%) are lower than for 
C IV. However, we cannot immediately conclude that the 
global covering factor of C iv absorbers is higher than 
those of N V and Si iv, because all of these values are just 
lower limits since some intrinsic NALs could have full cov- 
erage and also because we were not able to search for 
N V NALs over a wide range in VgMft- If we consider all 
three transitions, 43% (16 out of 37) of the quasars have 
at least one class A intrinsic NAL (in at least one of the 
three transitions) and 54% (20 out of 37) of the quasars 
have at least one class A or class B intrinsic NAL. 

4.5. Distribution of Coverage Fractions 

The normalized distribution of the coverage fraction 
for individual NAL components is presented in Figure [HI 
where it is broken down by transition and by reliability 
class. We plot only components for which we get physical 
covering fractions. The mean and median values of C/ in 
class A+B NALs are, respectively, 0.38 and 0.33 for N v, 

0. 50 and 0.48 for C iv, and 0.56 and 0.52 for Si iv. Putting 
together NALs from all three transitions, we find mean 
and median values of Cf of 0.48 and 0.45, respectively. 
Focusing on the Cf distribution among C iv NALs (the 
largest of all populations shown in Fig. we note that in 
class A-l-B NALs are distributed uniformly between and 

1, while more than half of the class C NALs have coverage 
fraction larger than 0.9 (this was one of the criteria for 
placing a component in this class). 

In addition to the NAL components shown in Figure [9l 
Tabled] shows that 56%, 34%, and 58% of C iv, N v, and 
Si IV NAL components in the homogeneous sample have 
unphysical measured Cf values. This usually indicates 
a doublet ratio slightly greater than two. If we assume 
that all components whose Cf values are unphysical or 
are consistent with full coverage (i.e., Cf + cr{Cf) > 1.0) 
really have fuU coverage, we find that 73% (355/483), 41% 
(17/41), and 72% (131/182) of C iv, N v, and Si iv com- 
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ponents cover the background flux sources. If we con- 
sider only intrinsic systems, these fractions would be 59% 
(44/75)), 26% (8/31), and 31% (8/26), respectively Gen- 
erally, we display only components whose C/ values are 
physical and whose errors are small enough (i.e., tT(C/) < 
0.1) to avoid any uncertainties from Cj value estimation. 

For some intrinsic NAL systems, coverage fractions have 
been determined for more than one transition (e.g., Gan- 
guly et al. 2003; Yuan et al. 2002). G99 found that the 
C IV, N V, and Si iv coverage fractions are similar to each 
other for a few systems, while Petitjean & Srianand (1999) 
and Srianand & Petitjean (2000) noted that higher ioniza- 
tion transitions tend to have larger coverage fraction than 
lower ionization transitions. This type of effect could oc- 
cur if the sizes of the background sources (i.e., continuum 
source, BELR) and/or the absorber are not the same for 
all transitions. 

Unfortunately, our results do not allow a direct compar- 
ison between different transitions within the same system 
because (a) detecting multiple transitions from the same 
system was rather rare and (b) multiple components in the 
same intrinsic NAL often have different C/ values, pre- 
venting us from assigning a single C/ value to a NAL or 
a system. Therefore, we compare the overall distributions 
of the physically meaningful (between and 1) C/ values 
of the 34 C IV, 25 N v, and 20 Si iv components found in 
23 C IV, 9 N V, and 9 Si iv class A and B intrinsic NALs, 
in Figure [51 An alternative illustration of these distribu- 
tions is shown in Figure [TO] where we plot the coverage 
fraction of each NAL component against its offset veloc- 
ity. A visual inspection of these two figures suggests that 
the associated intrinsic N v NALs prefer smaller values of 
Cf than the C iv and Si iv intrinsic NALs. Nearly all N v 
components show coverage fractions less than 0.5, while 
the coverage fractions of C iv and Si iv NALs range from 
nearly up to nearly 1. This difference is also reflected in 
the mean and median values of C/ reported above. How- 
ever, a Kolmogorov-Smirnov test yields a probability of 
17% that the C iv and N v C/ distributions were drawn 
from the same parent population. Similarly, the probabil- 
ity that the N v and Si iv Cf distributions were drawn 
from the same parent population is 7%. Larger samples 
are needed to reach a definitive conclusion. 

4.6. Other Relations between Intrinsic NAL and Quasar 

Properties 

With our large sample of intrinsic NALs we can search 
for relations between the NAL properties (e.g., Cf vs 
shift and rest equivalent width vs v shift) and the prop- 
erties of the quasars hosting them {Lopt, Lj-adio, TV). Fig- 
ure [TT] shows plots of these parameters against each other, 
including class A-l-B NALs. For some quasars, we plot 
just the upper limits on Lradio, since the corresponding 
radio sources have not been detected. We do not find any 
significant correlations in these plots. We also searched 
for a correlation between the offset velocity and the cover- 
age fractions and rest equivalent width of C iv, N v Si iv 
NALs but we did not find any (see Figure [TO)) 

4.7. Ionization State 

Traditionally, BALs have been classified according 
to their ionization level into HiBALs (high-ionization 



BALs), LoBALs (low-ionization BALs), and FeLoBALs 
(extremely low-ionization BALs that show Fe ii lines; e.g., 
Weymann et al. 1991). Motivated by the classification of 
BALs and by the work of Bergeron et al. (1994) on clas- 
sifying intervening NALs, we attempt to classify intrinsic 
NALs (classes A plus B) in a similar manner. We consider 
the range of conditions in the systems presented in this 
paper in conjunction with our concurrent study of lower 
redshift intrinsic NALs in the iTST/STIS Echelle archive 
(Ganguly et al. 2007; in preparation). This low redshift 
sample adds an important element because the O vi dou- 
blet can often be examined, unlike in the spectra of z = 2- 
4 quasars of this paper for which the Lya forest is quite 
thick. 

We find that we can classify the intrinsic NALs in our 
sample into 2 categories according to the relative absorp- 
tion strengths in various transitions: 

Strong C iv NALs. — In these systems Lya is often sat- 
urated and black. N v is detected in some of these 
C IV systems, but it is weaker than C iv. 

Strong N v NALs. — The corresponding Lya lines have 
an equivalent width that is less than twice the 
N V A1239 equivalent width. C iv is detected in 
some of these N v systems, but is typically weaker 
than N v. 

To clarify, in the case of a NAL system which has both 
C IV and N v transitions, we classify it as a strong C iv 
system if Lya in the system is black or if the total equiv- 
alent width of Lya is larger than twice the N v equivalent 
width. Otherwise, we classify them as a strong N v sys- 
tem, even if C iv is stronger than N v. 

One possible explanation for the physical difference be- 
tween C IV and N v systems is a difference in ionization 
parameter, however the column density, the metallicity, 
and/or the metal abundance pattern could also affect these 
conditions. Photoionization models will be required for the 
conclusive discussion (Wu et al. 2007; in preparation). 

We applied the above criteria to 28 class A and 1 1 class 
B NALs, and classified them into 28 strong-C iv systems 
(17 class A and 11 class B) and 11 strong-N v systems (all 
of them are class A). The Zabs — 2.7164 system toward 
IIS1700-I-1055 has an O vi doublet which is much stronger 
than the N v, C iv, and Lya lines. Although we classi- 
fied it as an N V intrinsic NAL, it would likely fall in the 
additional category of a strong O vi system, if we allowed 
for that category. Finally, we added a subclassification 
based on the detection of low-ionization (IP < 25 eV) or 
intermediate-ionization (IP = 35-50 eV) absorption lines 
in the same system. We note that BAL systems some- 
times have extremely low-ionization transitions such as 
Mg II (IP = 15 eV; e.g., HaU et al. 2002). The results 
are summarized in Table [Jj where we report the transitions 
detected in each of the above systems. Notes on individual 
systems and velocity plots are included in the Appendix. 

Table [H shows that among the 28 C iv intrinsic systems, 
15 have low-ionization lines, while 25 have intermediate- 
ionization lines. Two systems (at Zabs — 2.7699 in 
Q1425-I-6039 and at Zabs = 2.3198 in Q1548-h0917) are 
sub-DLA systems, which casts some doubt on their intrin- 
sic classification (see the Appendix for a discussion of these 
individual cases). Nonetheless, we see that more than half 



10 



Misawa et al. 



of the intrinsic C iv NALs have related low ionization ab- 
sorption lines. Among the 11 N v systems, 5 have interme- 
diate ionization absorption, and only 1 has low ionization 
absorption. This is another noteworthy difference between 
C IV and N v intrinsic NALs. 

In 41% of intrinsic systems (16 out of 39) we detect 
low-ionization transitions, and 55% of quasars (11 out of 
20) that have at least one intrinsic system, have at least 
one system with low-ionization transitions. This is much 
higher than the fraction of LoBAL quasars among BAL 
quasars (17% according to Sprayberry & Foltz 1992, 13% 
according to Reichard et al. 2003). 

Unfortunately, O vi is in the Lya forest or is not cov- 
ered in our spectra, making it quite difficult to establish 
the ionization properties of our systems. However, an ad- 
ditional class of O VI associated NALs (possibly intrinsic) 
has emerged from our study of lower redshift quasars using 
the STIS/Echelle archive. These O vi systems are some- 
times redward of the emission redshift of the quasar, they 
have relatively weak Lya lines compared to O vi, and they 
do not always have partial coverage (Ganguly et al. 2006; 
in preparation). A few of our N v intrinsic NALs with 
non-black Lya may have very strong O vi absorption, and 
thus could fall in this third category. 

5. DISCUSSION 

5.1. Evolution of the Frequency of Intrinsic NALs With 

Redshift 

If we associate intrinsic NALs with quasar outflows, 
their evolution with redshift traces the evolution of quasar 
outflows. Moreover, since outflow properties depend on 
the properties of the accretion flow (at least in disk wind 
models), the redshift evolution of NALs provides con- 
straints on the evolution of quasar fueling rates. 

We begin by considering quasars with strong C iv AALs 
(Wrest > 1.5 A), which have a high probability of be- 
ing intrinsic because of their small velocity offsets from 
the quasars. At intermediate redshifts (z = 1-2), F86 
found strong C iv AALs in 33% of quasars, while GOl de- 
tected such a strong C iv AAL in only 2% of low-redshift 
quasars (z < 1), and suggested that this fraction increases 
with redshift. In our high redshift sample (z = 2-4), we 
only detect one rather strong AAL with Wrest ~ 1-0 A. 
This could mean that the frequency of strong AALs peaks 
at z ^ 2, which coincides with the redshift at which the 
quasar density peaks. There are, however, a number of 
noteworthy caveats to this apparent evolutionary trend, 
namely (a) the samples used at different redshifts have dif- 
ferent optical and radio luminosity distributions (the GOl 
sample includes a significant fraction of low-luminosity 
objects, for example), (b) strong AALs appear to prefer 
steep-spectrum radio-loud quasars (F86; Anderson et al. 
1987) which make up different fractions of the above sam- 
ples (the F86 sample is predominantly radio-loud), and 
(c) strong AALs need not always be intrinsic, in particu- 
lar they could come from galaxies in a cluster surrounding 
the quasar, which are more often found around radio-loud 
quasars (e.g., Yee & Green 1987, Elhson, Yee, & Green 
1991; Yee & Ellison 1993; Wold et al. 2000). 

We now turn our attention to the evolution of the frac- 
tion of quasars with intrinsic AALs, including weak ones, 
down to a rest-frame equivalent width of ~ 50 mA (us- 



ing samples constructed on the basis of partial coverage or 
variabihty). At z < 1.5, Wise et al. (2004) found that a 
minimum of 27% of quasars have intrinsic AALs, while at 
z ~ 2, Narayanan et al. (2004) found that a minimum of 
25% of quasars host intrinsic AALs. Similarly, G99 find 
that at least 50% of z ~ 2 quasars have intrinsic AALs, and 
in this study, we find that a minimum of 23% of quasars at 
z — 2-4 have intrinsic AALs. These fractions are consis- 
tent with each other, especially considering that some of 
the samples are rather small (only six quasars in the G99 
sample) . An important caveat here is that all of these re- 
sults refer to lower limits since NALs that do not show 
variability or partial coverage can still be intrinsic. This 
caveat not withstanding, the fraction of quasars hosting 
intrinsic AALs does not show significant redshift evolu- 
tion. A closely related question is the redshift evolution 
of the fraction of AALs that are intrinsic. At z < 1.5 this 
fraction is > 21% (Wise et al. 2004), at z ~ 2 it is > 23% 
(Narayanan et al. 2004) or > 60% (G99), and at z = 2-4 
it is > 27% (this work). Thus, this quantity does not ap- 
pear to evolve with redshift either, but this conclusion is 
also subject to the caveat noted above. 

5.2. Properties of Intrinsic Absorbers in Our Sample 
5.2.1. Range of Coverage Fractions 

In non-associated regions (i.e., |ws/ii/t| > 5,000 km s^^), 
we found that 10-17% of C iv NALs are intrinsic to the 
quasars. Richards et al. (1999) and Richards (2001) sug- 
gested that this fraction is as large as 36%, based on the 
excess of absorbers at > 5000 km s""'^ in steep spectrum as 
compared to flat spectrum radio- loud quasars. Our esti- 
mate is strictly only a lower limit, because not all intrinsic 
NALs can be identified using partial coverage. If we take 
both our result that 10-17% of C iv NALs are intrinsic and 
the Richards et al. number of 36% at face value and try 
to reconcile them, we conclude that only 30-50% of intrin- 
sic C IV NALs exhibit detectable partial coverage in our 
sample. We caution that the two samples were selected in 
different ways, leading to biases. For example, the sample 
of Richards et al. (1999) has a larger fraction 50%) 
of radio-loud quasars than our sample (32%). However, 
since they found that intrinsic NALs prefer radio-quiet 
quasars, the difference between our estimate of the frac- 
tion of quasars hosting intrinsic NALs and theirs would be 
larger than what it appears, leading to a larger number of 
fully covered intrinsic NALs. 

This conclusion is bolstered by the fact that the dis- 
tribution of measured coverage fractions of intrinsic C iv 
NALs spans the range from to nearly 1 (within uncer- 
tainties). The covering factors are presumably related to 
the physical sizes of the absorbing structures, which could 
span a significant range. In this sense there is nothing 
special about the size corresponding to a covering factor 
of 1. If there is a cutoff in size, there is no reason it would 
occur particularly at this value. Thus we strongly expect 
that many intrinsic absorbers with unit coverage fraction 
exist. 

5.2.2. Ionization Conditions 

The wide range of absorber ionization conditions pre- 
sented in (e.g., "strong N v" vs "strong C iv" ab- 
sorbers) may be a result of a wide distribution of distances 
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of the gas parcels from the continuum source and/or a 
wide distribution of their densities. The velocity offset of 
a system can serve as an indicator of its proximity to the 
continuum source under the assumption that parcels of gas 
in an outflowing wind arc accelerated outwards from the 
continuum source. Detailed photoionization modeling of 
selected systems can help us determine whether there is 
a relation between their ionization state of NAL systems 
and their apparent outflow velocity, which will allow us to 
determine whether there is ionization stratification in the 
outflow. The same models may also provide constraints 
on the density of the gas. In fact, photoionization mod- 
els for the associated absorber in QSO J2233— 606 (Gabel, 
Arav & Kim 2006) show that the higher-ionization kine- 
matic components are at lower blueshifted velocities than 
lower-ionization components. Moreover, preliminary re- 
sults of photoionization models for three of the quasars 
in our sample (Wu ct al. 2007, in preparation) suggest 
that absorbers at different velocities in the same quasar 
can have a wide range of densities. An alternative way 
of constraining the density is through monitoring obser- 
vations, which can set limits on the recombination time 
scale of the absorber, hence constrain the density (see, for 
example, the discussion in Wise et al. 2004, Narayanan et 
al. 2004, and M05). 

It is also important to ask whether the low ionization 
transitions in intrinsic C iv NALs are aligned with the 
high ionization transitions. An inspection of the velocity 
plots of individual systems in the Appendix shows that 
the line centers and profiles of low-ionization lines in in- 
trinsic NAL systems are almost always similar to those of 
high-ionization lines. This leads us to conclude that the 
low-ionization lines arise in the intrinsic absorbers and not 
in the ISM of the host galaxy. 

5.2.3. Intrinsic NALs at Large Velocity Offsets 

Our finding (supported by the independent method of 
Richards 2001) that a significant fraction of NALs, even at 
high velocity offsets from the quasar redshift, are intrinsic 
(10-17% of C IV NALs and 15-20% of Si iv NALs) has im- 
plications for cosmology since it affects the estimates of the 
density of intervening systems per unit redshift [dN/dz). 
For example, at z = 2-4, dN/dz is estimated to be 3.11 for 
C IV and 1.36 for Si iv systems with a rest frame equivalent 
width greater than 0.15A (Misawa et al. 2002). If these 
values are reduced by 15 30% after correction for intrin- 
sic systems, this would require some adjustment of models 
for cosmological evolution of these systems. Furthermore, 
the contamination of the Lya forest by very high metallic- 
ity intrinsic NALs could seriously bias estimations of the 
metallicity of the forest. 

5.3. Implications for Models of Quasar Outflows 

Th(^ rc^sults of our survey have direct implications for 
models of quasar outflows, whatever the acceleration 
mechanism. 

1. The minimum fraction of quasars with at least 
one intrinsic NAL (C iv, N v, or Si iv) is 43- 
54%. If there is an additional family of intrinsic 
O VI absorbers, without corresponding C iv or N v 



lines, the minimum fraction of quasars with intrin- 
sic NALs will increase further. This leads to a con- 
straint on the solid angle subtended by the NAL 
absorbers to the background source(s). We can in- 
terpret this result in the context of a simple (but 
clearly not unique) picture in which all quasars in 
our sample are similar to each other and the NAL 
gas is located at intermediate latitudes above the 
streamlines of a fast, equatorial outflow (a BAL 
wind). This geometry is similar to that depicted 
in Figure 13 of GOl, although it needs not apply 
specifically to an outflow launched from the accre- 
tion disk; a pressure-driven outflow launched from a 
distance of an order of 1 pc from the black hole may 
have a similar geometry. In this picture, the above 
observational constraint translates directly into the 
opening angle of the NAL zone, given the opening 
angle of the BAL wind. For an opening angle of 
the BAL wind of approximately 7-12° (10-20% of 
quasars are BAL quasars; e.g., Hamann, Korista, 
& Morris 1993) the angular width of the NAL zone 
turns out to be 25-27°. In this picture the polar re- 
gion, at latitudes above the NAL zone may be filled 
with highly- ionized gas. This idea is consistent with 
the finding of Barthel, Tytler & Vestergaard (1997) 
that strong associated C iv NALs with Wrest > 3 
A are most common at angles far from the jet axis 
in radio-loud QSOs (as high as 45°). Such strong 
NALs may represent intermediate viewing angles 
between intrinsic NALs and BALs. Other simple 
interpretations are also possible: for example, the 
NAL gas may be distributed isotropically around 
the central engine, or some quasars may be more 
likely than others to host such absorbers (e.g., the 
width of the NAL zone may depend on luminosity, 
as suggested by Elvis 2000). 

2. The density of intrinsic (class A-|-B) C iv NALs 
per unit velocity interval, dN/df3, increases from 
4.1 at Vshift < —5,000 km s~^ to 18 at Vshift > 
—5,000 km s~^. This is a significant increase with 
a chance probability of only 5%. These results 
have a number of implications for outflow models. 
First, the high-velocity NALs (71% of class A-hB 
NALs at Vshift < —5,000 km s~^) are unlikely to 
be associated with pressure-driven outflows because 
these outflows can only reach velocities of order 
10^ km s"^ The distribution oidN/dp with /3 does 
not provide a straightforward constraint on the ac- 
celeration mechanism because the models do not 
predict how individual parcels of NAL gas are dis- 
tributed. One could take the higher value of dN/d(3 
at low velocities as an indication of the existence 
of a separate population of absorbers, which could 
perhaps be identified with a pressure driven wind. 
On the other hand, the velocity of a magnetocen- 
trifugal or line-driven accretion disk wind increases 
smoothly with radius, which may lead to higher 
density of NALs per unit velocity at low velocities 
(assuming that the gas parcels are condensations 
that form near the base of the wind and are accel- 
erated outwards). Economy of means and Occam's 
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razor favor the latter explanation since it applies 
to both high- and low-velocity NALs and invokes 
a single acceleration mechanism. This issue can be 
addressed observationally by constraining the dis- 
tance of the low- velocity (associated) NAL gas from 
the central engine. This can be done by combining 
constraints on its ionization parameter (obtained 
from photionization models) with constraints on its 
density (obtained from variability studies). 

3. The C IV column densities derived from fitting the 
line profiles, can yield a preliminary estimate of the 
total hydrogen column densities (H i-l-H ii) in the 
absorbers. More specifically, we can apply the cal- 
culations of Hamann (1997) to estimate an upper 
limit to the hydrogen column density of absorbers in 
our sample from our highest measured C iv column 
density of log Nc iv — 15.86 (in the Zats = 2.5597 
NAL in HE0130-4021). If we assume an abundance 
of C IV relative to C of log [Nciv/Nc] w -0.34 
(this limit corresponds to the optimal ionization 
parameter that maximizes the C iv abundance; see 
Figure 2 of Hamann 1997) and a solar abundance of 
C relative to H, log (C/H) = -3.44 (Grevese & An- 
ders 1989), we obtain A^h ~ 4 x lO^^ cm"2. if the C 
abundance is super-solar (e.g., Hamann 1997), this 
limit will become smaller, of course. Higher total 
column densities have also been observed in NALs; 
for example, Arav et al. (2001a) find log A^h > 20.3 
in QSO J2359-1241. Moreover, total column den- 
sities in BALs can be an order of magnitude higher 
(e.g., in PG 0946+301; Arav et al. 2001b). 

We emphasize that the column densities estimated 
above refer only to the gas responsible for the C iv 
NALs. This gas could well have the form of fila- 
ments embedded in a hotter medium, as predicted 
by many models (see fjl]and references therein). If 
that is the case, then the total column density of 
the outflow can be considerably larger than what we 
have estimated above. In fact. X-ray spectroscopy 
of nearby AGNs does indicate a multi-phase ab- 
sorber structure in which the column density of the 
hot medium is higher than that of the cold medium 
by a factor of 30 or more, (e.g., Netzer et al. 2003, 
Kaspi et al. 2004, and references therein) 

For the specific case of radiatively-accelerated outflows, 
we can ask whether the observed high velocities of some 
NALs are attainable in the context of the model. Hamann 
(1998) derives the following expression for the terminal 
velocity of a radiatively accelerated wind 

= 9, 300 r-y^ (^^^ - 0.1 Ms) ' km s'^ , (13) 
V JV22 / 

where r = 0.1 ro.i pc is the distance of the absorber from 
the continuum source, L — 10^^ L46 erg s~^ is the quasar 
luminosity, Mbh ~ 10^ Afg M0 is the black hole mass, 
and / = 0.1 /o.i is the fraction of continuum photons ab- 
sorbed or scattered by the gas, and Nh = 10^^ A22 cm^^ 
is the column density of the absorbing gas. Noting that 
Li^/M^ « 0.77 (L/LEdd) (where Lsdd is the Eddington 
luminosity) and that ro.i/Mg « 2 ^4 (^4 = Vg and 



rg = GMbh / is the gravitational radius), we may recast 
equation p3| as 

Voo = 5, 000 ^-^''^ - 0.13) ' km s-i . (14) 

\iV22 LEdd J 

As Hamann (1998) argues, /o.i ~ 1 for BAL quasars. The 
radiation-pressure dominated part of the accretion disk 
extends up to ^4 « 0.25 for parameters scaled as above 
(see Shakura & Sunyaev 1973). According to most mod- 
els for accretion disk winds (Murray et al 1995; Proga et 
al. 2000; Everett 2005), the inner launch radius of the 
wind corresponds to ^4 « 0.06-0.25. Moreover, the C iv 
column densities that we find observationally (see discus- 
sion above) indicate total hydrogen column densities in 
NALs corresponding to A''22 ~ 0.01. Under these condi- 
tions, radiation pressure can accelerate the NAL gas to 
terminal speeds well in excess of 10^ km s^^, and the ob- 
served blueshifts of intrinsic NALs can be easily explained. 

There is an important caveat to the above estimate. We 
have assumed that the column density that we have mea- 
sured from the C iv NAL profiles traces all of the matter 
that is being accelerated. But, as we note earlier in this 
section, the NAL gas may be embedded in a hot medium 
and it may represent only a small fraction of the total 
gas mass. In such a case, the attainable terminal velocity 
will decrease accordingly. More specifically, if the NAL 
gas represents up to 10% of the total column density, then 
equation (fT4)l can still reproduce the observed NAL veloci- 
ties. But if the total column density is more than an order 
of magnitude higher than the C iv NAL column density, 
then a different acceleration mechanism must be sought. 

6. SUMMARY AND CONCLUSIONS 

We have constructed a large, relatively unbiased, equiv- 
alent width limited sample of intrinsic narrow absorp- 
tion line (NAL) systems found in the spectra of z = 2-4 
quasars. This sample comprises 124 C iv, 12 N v, and 
50 Si IV doublets, which were separated from interven- 
ing NALs on the basis of their partial coverage signature. 
After assessing the reliability of the determination of the 
intrinsic nature of these NAL systems, 28 are deemed reli- 
ably intrinsic ("class A"), 11 are deemed possibly intrinsic 
("class B") and 111 are deemed intervening (or unreliable 
intrinsic candidates; "class C"). Using this sample of NAL 
systems, we study their demographics, the distribution of 
their physical properties, and any relations between them. 
We also consider the implications of these results for mod- 
els of outflowing winds. Our findings and conclusions are 
as follows: 

1. The fraction of intrinsic C iv systems in our sam- 
ple NALs is 11-19%. This value increases to 33% 
if only associated systems, within 5,000 km s~^ of 
the quasar emission redshift, are considered. This is 
roughly consistent with previous studies of associ- 
ated C IV NALs at different redshifts and employing 
different methods. It is important to note, however, 
that all of the above fractions arc, in fact lower lim- 
its to the true fractions. The fraction of intrinsic 
Si IV systems is 14-18%, although all of these are 
found in the non-associated regions of the spectra, 
at offset velocities VsMft < —5,000 km s^^ relative 
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to the quasar redshifts. Severe contamination by 
the Lya forest allowed us to search only for associ- 
ated N V systems. We found 75% of such systems 
to be reliably intrinsic. 

2. The minimum fraction of quasars that have one or 
more intrinsic NALs is 43-54% . While C iv intrin- 
sic NALs were detected in 24-32% of the quasars 
with adequate velocity coverage, N v intrinsic NALs 
were detected in 19% of the quasars, and Si iv in- 
trinsic NALs in 19 24% of the quasars. These are 
lower limits because our spectra do not have full 
offset velocity coverage for these transitions and be- 
cause some intrinsic absorbers may not exhibit the 
signature of partial coverage. This places a con- 
straint on the solid angle subtended by the absorber 
to the background source(s). 

3. We find that 10-17% of non-associated C iv NAL 
systems are intrinsic. In cosmological applications, 
non-associated systems are typically taken to be in- 
tervening. Thus, our result is important because it 
shows that it is necessary to correct derived cos- 
mological quantities for a contamination by intrin- 
sic NALs. A similar conclusion was reached by 
Richards et al. (1999) and Richards (2001) based on 
a statistical study of NALs in different quasar sam- 
ples, although they estimated a somewhat higher 
contamination of ~36%. Taking the two estimates 
at face value (i.e., neglecting the possibility that 
they differ because of systematic effects and de- 
manding that they should be reconciled) we are led 
to the conclusion that only 30-50% of intrinsic C iv 
NALs exhibit partial coverage. 

4. The coverage fractions of intrinsic NALs in our sam- 
ple, span almost the entire range from to 1. Since 
there is a range of sizes for intrinsic gas parcels, this 
also implies that there are a significant number of 
intrinsic absorbers that have full coverage, and thus 
were not detected in our survey. There is no ap- 
parent relationship between coverage fraction and 
velocity offset. Nor is there any relation between 
the NAL properties or frequency of incidence and 
the properties of the host quasars (such as optical 
or radio luminosity). 

5. We consider the ionization structure of the 39 class 
A and B intrinsic NAL systems in our sample, and 
find two major categories, which may represent ab- 
sorbers of different densities and/or at different dis- 
tances from the source of the ionizing continuum. 

(a) "Strong C iv" systems are characterized by 
strong, partially covered C iv doublets, strong, usu- 
ally "black" , Lyo! lines, and relatively weak or un- 
detectable N V doublets. Of the 28 systems in this 
category, 25 also have intermediate ionization lines 
detected (such as Si iii , C iii , and/or Si iv ) and 
15 systems also have low ionization absorption de- 
tected (e.g., O I, Si II, Al ii). In cases where the Lya 



profile is black, it is clear that it does not arise in 
the same gas parcel as the C iv absorption, since 
the coverage fraction for Lya is clearly 1. These 
C IV systems cannot be distinguished from typical 
intervening systems except by partial coverage. 

(b) "Strong N v" systems are characterized by 
strong N v lines, and relatively weak, non-black 
Lya lines (with less than twice the equivalent width 
of their N v lines). C iv and O vi lines may also be 
detected in these systems, and in some cases O vi 
may be stronger than N v. We find 11 systems in 
this way, 5 of which have intermediate ionization 
transitions detected, and only 1 of which has low 
ionization transitions detected. 

6. About 53% of class A or B NAL systems include 
low ionization lines. This fraction is much higher 
than that of LoBAL quasars among all B AL quasars 
(13-17%). In the 15 C iv systems with detected low 
ionization lines, the line profiles of the C iv and low 
ionization lines are similar. In particular, low ion- 
ization lines are rarely detected at velocities other 
than those of the partially covered C iv compo- 
nents, implying that both families of lines arise in 
the same parcels of gas. 

7. Our detection of a significant population of intrin- 
sic, high- velocity NALs {vshift ^ lO^km s"^) dis- 
favors scenarios in which the absorbing gas is as- 
sociated with a pressure driven wind that does not 
originate very deep in the potential well of the black 
hole. The low-velocity NALs that we have detected 
could still originate in pressure driven winds. How- 
ever, economy of means leads us to prefer accretion 
disk winds because these can explain both the high- 
and the low- velocity NALs. 

Repeated observations of the NALs in this sample at 
the same spectral resolution and S/N would be particu- 
larly valuable for the following reasons. First, they will 
allow us to search for NAL variability, which can be used 
to confirm intrinsic NALs and to probe the relation be- 
tween variability and partial coverage of intrinsic NALs. 
Second, variability sampled via multi-epoch observations 
carries additional information that may allow us to con- 
strain properties of the absorber (e.g., density and distance 
from the continuum source; see applications by Hamann 
et al. 1997a; Narayanan et al. 2004; M05). Multi-epoch 
observations of such a large sample of NALs have never 
been carried out, in spite of their utility in providing cru- 
cial constraints on quasar outflows. The sample of NALs 
presented here is an ideal starting point for such work. 
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Sample Quasars and Their Properties 
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^Quasar names are based on B1950 coordinates. 
'^V magnitude from Veron-Cetty & Veron (2003). 

"^R magnitude from USNO-A2.0 Catalog (Monet et al. 1998), except for Q1330+0108, whose R magnitude comes from the USNO-B Catalog (Monet ct al. 2003). 

■^Optical flux density at 4400 A in the quasar rest frame (converted from observed flu 
Q1055+4611 were corrected for Lya forest contamination as described in ^of the text. 

''Observed radio flux density from NVSS (1.4GHz; Condon et al. 1998), Griffith et al. (1995; 4.85GHz), or Becker, White, & Edwards (1991; 4.85GHz). If 
no radio source is detected within lO" of the optical source, we use the detection limit of the survey as an upper limit to the radio flux. The radio fluxes of 
HE0130-4021 and Q1759+7539 come from Smith U Wright (1980) and Hook ct al. (1996), respectively. 

^Radio frequency corresponding to the observed radio flux density. 

^Radio flux density at 5 GHz in the quasar rest frame, obtained from the observed flux, assuming that/j^ oc v^ "^ . 
"^Radio loudness parameter, defined in ^of the text, 
^radio-loud or radio-quiet quasar (see ^for definition). 

■'Red limit of the velocity window in which we searched for C IV, N V, or Si IV NALs. 
^This lensed quasar is amplified by a factor of ~3.1 (Barvainis Ivison 2002). 
^This lensed quasar is amplified by a factor of 15.38 (Kormann et al. 1994). 
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Table 2 

Census of Poisson Systems, NALs, and Kinematic Components 
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'^Combined multiple NALs that lie within 200 km s ^ of each other. 
See discussion in fJ3]of the text. 

^Narrow kinematic components deblended by Voigt profile fitting with 
MINFIT. Numbers in parenthesis are number of components whose C/ 
values are physical (i.e., < C/ and C/ < 1). 

'^AU lines that are detected at greater than a 5cr confidence level, 
Wobs/criyVobs) > 5, regardless of their equivalent width. 

'^Doublets with rest-frame equivalent widths of their stronger (bluer) 
components larger than the minimum equivalent width that correspond 
to the 5cr detection limit, Wmin- These limits are Wmin (C iv) — 
0.056 A, W^in (N v) = 0.038 A, and iy„,„ (Si iv) = 0.054 A (see 
g3]of the text). 
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Table 3 

Statistical Properties of Poisson Systems of NALs 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 








CIV 






N V 






Si IV 








AAL 


non-AAL 


total 


AAL 


non-AAL 


total 


AAL 


non-AAL 


total 






0.64 


17. 


18. 


2.3 


0.32 


2.6 


1.3 


12. 


14. 






0.17 


4.9 


5.1 


0.58 


0.079 


0.66 


0.32 


3.3 


3.6 






3(3) 


11 (8) 


14 (9) 


9 (7) 


0(0) 


9 (7) 


0(0) 


7(7) 


7(7) 


class A 


dN/dz 


4.7 


0.63 


0.78 


3.9 


0.0 


3.4 


0.0 


0.57 


0.5 




dN/dl3 


18. 


2.3 


2.8 


16. 


0.0 


14. 


0.0 


2.2 


2.0 






3(3) 


20 (11) 


23 (12) 


9 (7) 


0(0) 


9 (7) 


0(0) 


9 (9) 


9(9) 


class A+B 


dN/dz 


4.7 


1.2 


1.3 


3.9 


0.0 


3.4 


0.0 


0.73 


0.66 




dN/df3 


18. 


4.1 


4.6 


16. 


0.0 


14. 


0.0 


2.8 


2.5 






6(5) 


95 (34) 


101 (34) 


3 (3) 


0(0) 


3 (3) 


4(4) 


37 (24) 


41 (24) 


class C 


dN/dz 


9.4 


5.5 


5.6 


1.3 


0.0 


1.1 


3.2 


3.0 


3.0 




dN/d(3 


36. 


19. 


20. 


5.2 


0.0 


4.6 


13. 


11. 


11. 






9 (13) 


115 (37) 


124 (37) 


12 (37) 


0(0) 


12 (37) 


4 (24) 


46 (37) 


50 (37) 


Total 


dN/dz 


14. 


6.6 


6.9 


5.2 


0.0 


4.6 


3.2 


3.7 


3.7 




dN/d(3 


54. 


24. 


25. 


21. 


0.0 


18. 


13. 


14. 


14. 


NALs'^ 


A/Total 


33% 


10% 


11% 


75% 




75% 


0% 


15% 


14% 




A+B/Total 


33% 


17% 


19% 


75% 




75% 


0% 


20% 


18% 


Quasars'* 


A/Total 


23% 


22% 


24% 


19% 




19% 


0% 


19% 


19% 




A+B/Total 


23% 


30% 


32% 


19% 




19% 


0% 


24% 


24% 



"The total redshift and speed intervals considered in the determination of dN/dz and dN/df} 

^N denotes the number of Poisson systems (i.e., groups of NALs that lie within 200 km of each 
other; see ijSj. The number in parenthesis gives the number of quasars in wrhich these Poisson systems are 
found. 

'^Percentage of class A and class A-fB NALs relative to all NALs, broken down by transition and by 
velocity window (associated vs non-associated). 

'^Percentage of quasars hosting class A and class A-l-B NALs relative to all quasars, broken down by 
transition and by velocity window (associated vs non-associated). 
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Table 4 



Ionization Levels of Class A and Class B Systems" 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 






Low-Ion.'' 


Interm.-Ion.° 


Civ 


N V 


VI 


Ionization 


QSO 


^abs 


13-24 eV 


33-48 eV 


65 eV 


98 eV 


138 eV 




Class'' 


Class A Systems 


HFm 4.091 

X1.IZjUJ.OU — 4:UZ J. 




i\r 


V 




IN 




1 

J., 


TTT 




9 Q7A.Q 


N 






Y 


V 


9 


TI 
1 1 




^ 97S1 


N 


N 




V 




9 
z, 


TI 
1 1 




3.2818 




Y 




Y 




2 


HI 


on4'in— 1 sin 


2.2307 


N 


Y 




Y 




2 


HI 


V^UOoO-rOoUi 




N 

i\ 


V 


V 






1 

1 , 




00*=; /1 9-1- /I /I 


9 Q791 


N 
i\ 


V 


V 






1 

1 , 






9 A^AA 


IN 


V 


V 

I 


IN 




1 

1 , 


TTT 
ni 




9 fit^l 7 
Z.DOl / 


Y 
1 


V 

1 


V 

1 


V 




1 


TTTT 




9 SPiSQ 
Z.oOoy 


IN 


IN 


V 

1 


V 

1 


1 


9 


XJ 
XI 






V 


V 




N 




1 

± , 


TTTT 

XlXXj 






IV 


V 




V 

1 


V 


9 
z, 


TIT 
XIX 


Q1017+1055 


«2.97 


Y" 


Y" 


Y^ 






2, 


HIL 




fa3.03 


N 


N 


Y" 






2, 


H 




r;3.11 


N 


Y 




Y^ 




2, 


HI 


Q1055+4611 


3.5314 


Y 


Y 








1, 


HIL 


Q1107+4847 


2.7243 


Y 


Y 


Y 


N 




1, 


HIL 


Q1330+0108 


3.1148 


Y 


Y 


Y 






1, 


HIL 


Q1425+6039 


2.7699 


Y 


Y 


Y 






1, 


HIL (DLA) 


Q1548+0917 


2.3198 


Y 


Y 


Y 


N 




1, 


HIL (DLA) 




2.6659 


Y 


Y 


Y 


N 




1, 


HIL 




2.6998 


Y 


Y 


Y 


N 


N 


1, 


HIL 


Q1554+3749 


2.3777 


Y 


Y 




N 




1, 


HIL 


HS1700+6416 


2.7125 


N 


N 


Y 


Y 




2, 


H 




2.7164 


N 


N 


N 


Y 


Y 


2, 


H 


HS1946+7658 


2.8928 


Y 


Y 


Y 


N 




1, 


HIL 




3.0385 


N 


N 


Y 


N 




1, 


H 




3.0497 


Y 


Y 


Y 


Y 




1, 


HIL 


Class B Systems 


HE0130-4021 


2.2316 


N 


Y 


Y 






1, 


HI 


Q0241-0146 


3.0451 


N 


Y 


Y 






1, 


HI 


HE0940-1050 


2.8347 


N 


Y 


Y 






1, 


HI 


Q1017+1055 


2.5408 


N 


N 


Y 


N 




1, 


H 


Q1055+4611 


3.3658 


N 


N 


Y 






1, 


H 


HS1103+6416 


1.8874 


N 


Y 


Y 


N 




1, 


HI 




1.8919 


Y 


Y 


Y 


N 




1, 


HIL 


Q1334-0033 


2.2010 


N 


Y 


Y 


N 




1, 


HI 


Q1548+0917 


2.6082 


Y 


Y 


Y 


N 




1, 


HIL 


HS1700+6416 


2.4330 


Y 


Y 


Y 


N 




1, 


HIL 




2.4394 


N 


Y 


Y 


N 




1, 


HI 



"Y" indicates that the system contains one or more lines of a given type are detected. A "N" 
indicates that no lines from that category were detected, while a blank entry means that either the 
spectrum did not cover those transitions or the transitions were severely blending with Lya forest or 
data defects. 

''Low-ionization (from ions with ionization potentials between 13 and 24 eV): O I A1302, Si II A1190, 
A1193. A1260, A1527, Al ii A1671, C ii A1335. 

"^Intermediate-ionization lines(from ions with ionization potentials between 33 and 48 eV):Si III A1207, 
Si IV A1394, A1403, C ill A977. 

''l=C IV dominated system, 2=N v dominated system. H, I, and L mean that the system contains 
high-, intermediate-, and low-ionization lines, respectively. 

"These lines in Q1017-I-1055 are classified as "mini-BALs" because they have widths of 500- 
1,000 km s~^. See detailed discussion in the Appendix. 



Table 5 

Properties of Narrow Absorption Lines 



(1) 


(2) (3) 


(4) 


(5) 


(6) (7) (8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 










(A) 


log AT 


cr(v) 




Reliability 




Other 


QSO 


2e™ AAL" 


Ion 


(A) 


(cm-^) 


(km s"^) 




Class" 


Sample'^ 


lons^ 



Q0004+1711 2.8 



C IV 4843.0 2.7768 2.1281 -64375 

2.1271 -64473 

2.1275 -64431 

2.1280 -64385 

2.1283 -64357 

2.1286 -64334 

2.1289 -64309 

2.1290 -64298 



13.71±0.22 
14.21±0.05 
14.48±0.04 
12.44±1.80 
14.28±0.05 
12.74±1.32 11.0±13.1 
13.88±0.06 10.8±0.7 



82.0 

19.3±2.9 0.49^017 
11.9±0.7 0.97t°o°ol 
22.1±1.6 
7.4±13.4 
15.3±3.3 



CI 



L,H 



Lya(C II A1335, Si III A1207, 
Si IV) 



U-yo-0.07 
1.00 
1.00 
1.00 



1.00 



+ 0-08 
-0.07 



Si IV 4903.0 1.0584 2.5178 -30068 

2.5171 -30132 

2.5177 -30082 

2.5182 -30036 



42.6 

12.30±0.12 2.7±2.9 
13.97±0.07 11.5±0.5 
13.52±0.07 11.4±0.6 



1.00 
99+° ™ 

• +0.09 
-0.09 



CI L,H Lya(Si II AA1190,1193, 

Si II A1260, Si III A1207, 
Si IV) 



0.85: 



bO 

o 



Associated absorption line, \vshift\ ^ 5,000 km s ^. 

^Wavelength of flux-weighted line center. 
^Observed frame equivalent width of blue member. 
'^Redshift of flux weighted line center. 
^Velocity offset from quasar emission redshift. 
* Coverage fraction. 

^Partial coverage reliability classes, as defined in j|3.1l 

^Subsamples: L— radio-loud, Q— radio-quiet. A— associated (within 5,000 km of z^m), and H— homogeneous/complete sample with Wrest /(^iWrest ^ 5). 

^Other lines that are detected in the system. Lines in parenthesis are in Lya forest, i.e., less reliable. 

■'Evaluated with red member, because blue member blends with other lines, data defects, or echelle order gaps. 

^Lya, and N V seem to have disappeared ("normalized out"). 

^Voigt profile fit is not applied because of smooth wide absorption profile. 

"^Approximate value. 

"Fitting assuming C/ — 1. 
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Fig. 1. — Distributions of optical (4400 A) and radio (5 GHz) luminosities of the 37 quasars in our sample. The shaded bins in the lower 
panel correspond to upper limits to the radio flux, as listed in Table [T] Radio-loud objects cluster at 5GHz luminosities between 10^^ and 
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Fig. 2. — Two point correlation functions for C IV (solid line), N V (dashed line), and Si IV (dotted line) NALs in velocity offset space. 
Strong clustering is seen at Av <200 km s~^ (indicated by arrows in the main panel), especially for C IV NALs. Although we count only the 
number of pairs which are just next to each other, we get a similar result if we count all combinations. 
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Broad line profile? 




¥ 
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Line loc3kiiig? 



B2 



- model fit ia acceptable? 

- continutun fit is reliable? 

- without critical data defect? 





all Y 
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C3 





After rDmoving all unclassif iable ccunponents ! 
1* weak componenta at the edge of syatema 
2. weak components between stronger components 
3* heavily blended region 

(except for highest S/W regions} 
4, unphyaical componertts (Cf-1) 
No components left? 

N 



C2 



MINFIT gives Ci_+3af )<1 . 
for at leaet one component 



N 



A2 



MINFIT gives Cj.+tj(Cj^ )<1 . 
and 

pixel-by-pixel method gives C£+(J(Cf )*il . at the 



bottom of component 
for at least one component 





¥ 




N 
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CI 





Fig. 3. — Flow chart illustrating our scheme of NAL classification. The classes indicate how reliably we can determine whether the NALs 
are intrinsic or not, based on partial coverage analysis. A more detailed discussion of this scheme is given in ^3. II of the text. 
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Fig. 4. — Velocity offset distribution of C IV, N V, and Si IV systems. For each transition we plot separately the distribution of class A 
systems (reliable), class A+B systems (reliable+possible), and class C systems (intervening/unclassified). 
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Fig. 5. — Velocity offset distribution of all C IV systems compared to the distribution of class A+B (reliably and possibly intrinsic) C IV 
systems. 
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Fig. 6. — Distribution of rest-frame equivalent width of all C IV systems, compared to the distribution of class A+B (reliably and possibly 
intrinsic) C IV NALs. The vertical dashed line marks our detection limit of 0.056 A. 




Fig. 7.- Distribution of the flux-weighted line width (<j{v)\ calculated by eqn. 4) of intrinsic NALs (class A and B) and interven- 
ing/unclassified NALs (class C). The filled stars indicate the average cr(ii)-values for intrinsic NALs (53.4 km s~^) and intervening NALs 
(44.5 km s"^). 
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2 4 

Number of intrinsic NALs per quasar 

Fig. 8. — Distribution of the number of intrinsic systems per quasar. The top panel includes only class A (reliable) systems, while the lower 
panel includes class A+B (reliable+possible) systems. In each panel we show the distribution of C IV systems as a solid line, the distribution 
of N V systems as a dashed line and the distribution of Si IV systems as a dotted line. We note that the number of intrinsic systems per 
quasars in these histogram is strictly a lower limit (see discussion in ^of the text). 
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Fig. 9. — Distribution of coverage fractions of C IV, N V, and Si IV NAL components. We plot only components whose Cf values are 
physical (i.e., < Cf and Cf < 1) and evaluated with high reliability (i.e., a{Cf) j 0.1). For each transition we plot separately the distribution 
of class A NALs (reliable), class A+B NALs (reliable+possible), and class C NALs (intervening/unclassified). Even though some class-C 
NALs have C/ < 1. These values are consistent with 1 within their la errors. 
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Fig. 10. — Relationship between the velocity offset of intrinsic NALs (class A+B) from the quasars and their coverage fraction (top) and 
rest-frame equivalent width (bottom). In the top figure, we plot only components whose Cf values are physical (i.e., < Cf and Cf < 1) 
and evaluated with high reliability (i.e., cr{Cf) < 1.0). The upper panel in each case presents a wide velocity region (from —80,000 to 
10,000 km s~^). The lower panel focuses on the narrower velocity window from —7,000 to 2,000 km s~^. Open circles denote C IV NALs, 
filled circles denote N V NALs, and asterisks denote Si IV NALs. 
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Fig. 11. — Plots of NAL parameters (i.e., Cf and v^i^iff) against quasar properties (i.e., Li/(4400 A), L^(5 GHz), and TZ). Only intrinsic 
NAL components (class A+B), whose Cf values are physical (i.e., < Cf and Cf < 1) and evaluated with high reliability (i.e., criCf) < 0.1), 
are plotted. For quasars that were not detected in the radio, we plot upper limits on L^(5 GHz) and TZ. Open circles denote C IV NALs, 
filled circles denote N V NALs, and asterisks denote Si IV NALs. 
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Fig. 12. — (Sample) Results of partial coverage analysis applied to 206 C IV, N V, and Si IV doublets. Velocities are defined relative to 
the flux- weighted center of the system, as given in the first line of each block of Table [5] The first two panels of each set show the profiles of 
the blue and red members of a doublet on a common velocity scale, with the model proflle produced by MINFIT superposed as a dashed line. 
The positions of the kinematic components making up the model are marked with upward arrows in the bottom of each panel. The third 
panel of each set shows the two profiles together, along with the resulting coverage fractions and their la error bars. The coverage fraction 
determined by the pixel-by-pixel method is plotted in the form of black filled circles, while the coverage fraction determined by the fitting 
method is plotted in the form of green open circles. 
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Fig. 13. — (Sample) Velocity plots of various lines detected in the class B NAL system at Zaba = 2.2316 in HE0130— 4021. 



